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Abstract: An old result (astro-ph/9905303) by Jacobson implies that a black hole with 
Schwarzschild radius r s acquires scalar hair, Q oc r 2 ^, when the (canonically normalized) 
scalar field in question is slowly time-dependent far from the black hole, df(f> — [iM v with 
[ir s <C 1 time-independent. Such a time dependence could arise in scalar-tensor theories 
either from cosmological evolution, or due to the slow motion of the black hole within an 
asymptotic spatial gradient in the scalar field. Most remarkably, the amount of scalar hair 
so induced is independent of the strength with which the scalar couples to matter. We argue 
that Jacobson's Miracle Hair-Growth Formula® implies, in particular, that an orbiting pair of 
black holes can radiate dipole radiation, provided only that the two black holes have different 
masses. Quasar OJ 287, situated at redshift z ~ 0.306, has been argued to be a double black- 
hole binary system of this type, whose orbital decay recently has been indirectly measured 
and found to agree with the predictions of General Relativity to within 6%. We argue 
that the absence of observable scalar dipole radiation in this system yields the remarkable 
bound | /i| < (16days) _1 on the instantaneous time derivative at this redshift (as opposed 
to constraining an average field difference, A</>, over cosmological times), provided only that 
the scalar is light enough to be radiated — i.e. m < 10 -23 eV — independent of how the 
scalar couples to matter. This can also be interpreted as constraining (in a more model- 
dependent way) the binary's motion relative to any spatial variation of the scalar field within 
its immediate vicinity within its host galaxy. 
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1. Introduction 

Ten years ago Ted Jacobson discovered Jacobson's Miracle Hair Growth Formula 0, which 
shows how black holes can acquire scalar hair if they are embedded into regions in which the 
relevant scalar fields evolve slowly with time. We argue in this paper that this observation 
can significantly change the behaviour of black-hole binaries, and provide very robust new 
observational constraints on the existence of time-dependent scalar fields in the cosmological 
past. Unlike other recent efforts to constrain scalar-tensor theories using black hole physics 
[g], our arguments do not rely on black-hole super-radiance. 

Why consider light scalars? Light scalar fields provide one of the few known ways to 
consistently modify gravity over large distances, and one of the main ways in which candidate 
quantum theories of gravity can differ from one another is the spectrum of bosons they predict 
at very low energies. Such fields, if they exist and are sufficiently light, can mediate long-range 
forces that can observably compete with gravity. Because of this many of their properties 
have been carefully studied over the years H, || ||, 0, ^] . 

Furthermore, scalars commonly arise in fundamental theories, although it is more unusual 
for them to be light enough to mediate forces over macroscopic distances. Scalars are rarely 
this light because quantum corrections famously tend to give large contributions to scalar 
masses, even if they would have been massless at the purely classical level. But in some cir- 
cumstances symmetries can protect against masses, such as if the scalar is a pseudo-Goldstone 
boson || for a spontaneously broken approximate symmetry. This is what protects the mass 
of an axion 0, 01, or other axion-like fields fll^ l. Alternatively, light scalars can arise as 
part of the low-energy limit of an extra-dimensional model, with mass protection coming 
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from extra-dimensional symmetries like higher-dimensional general covariance |l3|, |j"4]| . In 
particular, extremely light scalars with masses of order the present-day Hubble scale natu- 
rally arise in this way in extra-dimensional proposals to address the cosmological constant 
problem [15, 16 1. 

Once a light scalar is present it is very unlikely to be time-independent, since a small 
potential can easily drive a light scalar to roll. So it is natural to examine how the time evo- 
lution of ambient scalar fields can affect the properties of black holes and other astronomical 
objects. It is in this context that Jacobson's observation of scalar hair growth by black holes 
is so interesting. 

We find in this paper that the scalar hair Jacobson predicts for black holes can have 
potentially measurable implications, by modifying the post-Newtonian response of orbiting 
black holes. In particular, we find it provides an unexpected source of scalar dipole radiation 
for an orbiting black hole, which would have been unavailable if the black hole had no hair 
(as is usually supposed). 

Specifically, for canonically normalized scalar fields that approach <p — > cf)^ + [it asymp- 
totically far from a pair of orbiting near-Newtonian black holes (with masses M A and M B ), 
we find the fraction of power radiated into scalar fields to be proportional to Aa 2 , where 
Aa = 4Gfj,(M A — M B ). We apply this expression to the quasar OJ287, which has been 
interpreted as a double black hole binary, leading to a bound | /i| < (16 days) -1 on the in- 
stantaneous time-variation of any sufficiently light (see below for precisely how light) scalar 
fields at redshift z ~ 0.3. 



Others [17] have also sought useful ways to apply Jacobson's Miracle Hair Growth For- 
mula, through its potential implications for black-hole accretion Ijlfifl . Rather than using (as we 
do) the scalar hair to argue for dipole scalar radiation, these authors argue that in quintessence 
models primordial black holes can accumulate significant amounts of the quintessence field, 
potentially leading to observable effects in gravitational radiation. 

This paper is organized as follows. In section 2.1, scalar-tensor gravity is introduced, and 
the no-hair theorem is explicitly illustrated for the static spherically-symmetric solutions of 
the vacuum field equations. In section the mass M and scalar charge Q of an astrophysical 
body are defined, in terms of both the behaviour of the fields asymptotically far from the 
body, and the functional derivatives of the point-particle action of the effective field theory in 
which the internal structure of the body has been integrated out. In section 2^3, Jacobson's 
time-dependent black hole solution with scalar hair (the Miracle Hair Growth Formula) is 
described. An analysis of the appropriate point-particle action is carried out, and it is found 
that this type of black hole couples to the scalar field differently than a star with the same 
scalar charge, implying that the standard formulas in the scalar-tensor literature for the 
Post-Keplerian orbital corrections need not apply to a binary system of two such black holes. 
Section 2.4 then argues that standard formulae do apply for the leading-order rate of scalar 
dipole radiation emitted by a such a binary system. In section ^ the result of this calculation 
is applied to the quasar OJ287. Finally, section | summarizes our results and prospects for 
future work. 
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2. Black holes: hair and dipole radiation 



This section contains the main formal results: describing black holes in scalar-tensor theories, 
with both static and slowly varying asymptotic scalar configurations. It also briefly summa- 
rizes the standard formulae for the dipole and quadrupole radiation rate from binary pairs in 
the near- Newtonian limit. 

2.1 Scalar- Tensor gravity and black-hole balding 

The simplest variants of scalar-tensor gravity involve a single scalar field, ip, and the space- 
time metric, g^ u , with action: S = S g + S ma _f Here S m at describes how the fields ip and g^ v 
couple to other matter kinds of matter and S g describes the mutual interactions of ip and 
the metric. Writing S g in a derivative expansion and keeping only up to two derivatives, the 
most general expression is 

S g = - J d 4 x^ {uM+gV [T^n^ + g^d^d^ +•••}, (2.1) 

where the form taken by the functions U(ip), J~(ip) and Q{ip) characterize a particular scalar- 
tensor theory. Here IZ^ denotes the Ricci tensor 1 constructed from the metric g^ v . Of 
particular interest are scalar-tensor theories which satisfy the weak equivalence principle. 
This can be ensured by requiring that matter couplings, S ma ,t, depend only on g^ v and not 
separately on ip. 

In practice the functions J- and Q need not be specified in detail because they can 
be dramatically simplified by performing appropriate field redefinitions of the form g^ v = 
A 2 {ip) g^y and ip = ip(<p), with the result 



S g = - — ^ J d 4 xv^ \gT (?V + 2 8^8^) + WttGV(p) + ••■ 



(2.2) 



where G is Newton's constant. This leaves V(<p) as the only important function governing 
the long-distance properties of the scalar-gravity couplings. For those theories for which 5 ma t 
depends only on g^ u = A 2 {(p)g fJLV and not separately on tp, all of the direct matter couplings 
of the field cp are described by the function A(ip) := ^4[-0((^)]. For such theories the two 
functions A(cp) and V((p) completely define the low-energy predictions. Brans-Dicke theories 
H comprise the widely studied special case where V((p) = and A((p) = e aip , for constant a. 

In what follows we require no assumptions about the scalar-matter couplings, but our 
interest is in situations where the scalar potential, V(ip) = A 4 (ip) U(ip)\^^y is negligible 
(making ip a minimally coupled massless scalar). The vacuum field equations obtained from 
( |2.2[) when V — are very simple 

Kfu, + 28^8^ = 0, (2.3) 

n<p = -L ^ ( ^=5 gr%<p) = o . (2.4) 



1 Conventions: we use metric signature ( — h H — h) and Weinberg's curvature conventions [[L9| (differing from 
MTW §J only by an overall sign for the Riemann tensor). Unless specified otherwise, units are chosen 
throughout with h = C = ks = 1. 
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Time-independent, spherically symmetric solutions to these equations are known |2l], ||, 
p^ l, and given by 

ds 2 = -f(r) dt 2 + + g(r) (d6 2 + sin 2 9 dc/) 2 ) and ip = ip(r) , (2.5) 
fir) V / 

where 



m = [i 

9{r) = r 2 ( 1 - - ) (2.6) 



r . 

l-Q 



r 

and tp(r) = <poc + | In ^1 

and the constants a and g must satisfy a 2 + g 2 = 1. 

The Schwarzschild black hole corresponds to the choice a = 1 and q = 0, in which case 
= ^oo is constant, and the constant I is revealed as the Schwarzschild radius: I = 2GM. 
In this context the no-hair theorems [^] state that a = 1 and q = is the only physically 
allowed black-hole solution. Within the family defined by eqs. ( |2.6| ), the problem with the 
more general solutions is that (p always diverges at r = i whenever g / 0. Since the field 
equation ( f2~3| ) shows that the curvature scalar, 1Z, also diverges there, r = t is an honest-to- 
God curvature singularity, rather than simply marking the coordinate singularity associated 
with an event horizon. 

2.2 Masses and scalar hair 

We next pause to define scalar hair more precisely, for use in later sections. One way to define 
the (ADM) mass [24], ||], M, and the 'scalar charge' || 2f| |26|| , Q, is through the large-r 



2 



asymptotic form of the metric and the scalar field 

1 2GM 1 G Q /n^X 

9tt = -H 1 and </> = ^ 1 . (2.7) 

r r 

In terms of these definitions we see that the integration constants I and q can be traded for 
M and Q, with 

2GM = at and GQ = y . (2.8) 

For later purposes it is worth elaborating a bit on the extent to which this definition of 
Q provides a useful definition of 'charge'. The main observation is that the scalar equation of 



motion, eq. (2^), can be read as a conservation law: d^J^ = 0, for J M := y/—g d^ip. Although 
the quantity Q is not the charge — i.e. the spatial integral of J 1 — for this conserved 
current, it is 'conserved' in the sense that eq. ( |2.4D implies that the combination \f—g g rr p>' 



2 Notice our conventions here for Q follow those of ref. [ p6[ , and so differ by a sign from those of refs. 
who define ip — ipoo + Gui/r + ■ ■ ■ . 
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is independent of r (with cp' := d r cp), for any static spherically symmetric configuration. Q 
as defined above is equivalent to this r-independent quantity because 

y/=ij g rr <p' = f g (J shxO = ^1 - £\ r 2 ip' sin 6 = y sin 9 = GQ sin 6 . (2.9) 

Our interest in the rest of this paper is in applications computing the radiation rate and 
post-Newtonian corrections to orbits, for which it is only the centre-of-mass motion of the 
black hole that is physically relevant. Furthermore, all objects are located far enough from 
one another to allow their weak mutual interaction to be described by the far-field regime. 
In this case it is more useful to define mass and charge in terms of the effective point-particle 
action of the black hole which governs the evolution of its centre of mass p7| , 25, 29, |3(i| ], 

S p := J dT£p[cp,g fJ , v ,d ,cp,R , /3j s,z,- ■ ■] , (2.10) 

where C p is a scalar whose dependence on the background metric and cp, and their derivatives 
depends on the physical properties of the object of interest (e.g. planet, star, neutron star, 
black hole). This dependence is obtained (see below) by matching the physical predictions 
of this effective theory to the predictions of a more microscopic description of the object 



in question g, ^9[. In eq. ( 2.10 ) C p is evaluated along the point-particle world- line, z^(t), 
denoted as T, with r representing proper time (arc- length) measured along T. (C p could also 
depend on other variables for objects with internal degrees of freedom, like spin.) Over-dots 
denote the derivative with respect to r, as in := dz^/dr. 

This action has two important uses. First, it can be used as a source in the field equations 
for cp and g^, and expressions for C p can be chosen to make the resulting solutions agree 
with the far-field solutions computed from an explicit description of the relevant object. (For 
example, for a star these solutions can be compared to the far-field solutions obtained by 
matching to interior solutions to the equations of hydrostatic equilibrium inside the star.) 

Second, its variation with respect to z^(t) can be used to determine the trajectory taken 
by the object moving within its local geometric environment, including any applied fields cp. 
For these latter applications it is useful to organize C p into an expansion in derivatives of the 
background fields, 

S p = J dr[m(<^) +h(<p)z> i z v d fl <p d v cp + k(cp)d, x yd»cp+ •••] , (2.11) 



since it is often only the first term, C p ~ m(cp) that is relevant for practical applications 
The coefficients like m, h and k take values whose order of magnitude can be determined 
on dimensional grounds in terms of the mass, M, size, L, and time-scale T of the source in 
question; typically with m ~ 0(M), and k ~ 0(MI?) and h ~ 0(MT 2 ) and so on. (Here 
the characteristic time-scale is T ~ L/v, where v denotes the typical speed of objects within 
the source; v ~ c = 1 for relativistic sources like neutron stars.) For practical applications 
to post-Newtonian orbits it is usually the non-relativistic limit of the point-particle action 
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evaluated in the field of the other body that is relevant, which we obtain by expanding in 
powers of v. (Although the systematics of this expansion are interesting in their own right 
|2^ , [30| we do not follow these complications here.) 

Our immediate interest is in constraining the function, C p (<p), that describes the black- 
hole source that gives rise to the explicit scalar-tensor solutions given by eqs. ( |2.6| ), above. 
This may be accomplished by solving the field equations in the presence of the point-particle 
action, 5(S g + S p ) = 0, for instance 

dp(y^ g ra v <p)+4nG(j2\ 5 3 [x-z(s)}=0, (2.12) 



and comparing the far-field result to the solutions of eqs. ( |2.6[ ). We do not use the derivative 
expansion, eq. ( |2.11 ), when doing so because this can break down when evaluating S p at 



self-field configurations sourced by the object itself, due to the divergence of these self-fields 
at the position of the source. 3 

Taking a static source, and integrating eq. ( 2.12j ) at a fixed time over a spherical region 



centered on the source position implies in particular 

= ^ j> WW ^Tg /V = f(r) g(r) <f'(r) = GQ , (2.13) 



-G 6Cp 

5(p 



where the far-left-hand term is evaluated at the source position (which generically diverges 
and so must be renormalized) . The Einstein equations give a similar expression involving the 
stress energy derived from C p . For example, for widely separated, weakly interacting objects 
moving within a field that asymptotes at large r to a static constant configuration, (p^, and 
assuming C p ~ m(ip), these considerations reproduce the standard results [|| 

M = m((/9oo) and Q = —m'^oo) . (2-14) 

Notice that use of the weak, far-field expansion does not imply any assumptions about 
the nature of the underlying source, and in particular does not exclude it being a relativistic 
star or a black hole Q. However, in the special case that the source is constructed purely from 
non-relativistic matter, and when this matter couples to ip only through the Jordan- frame 
metric, g^ u , we would additionally know that m(</?)dr = modf, where f is measured using 
the metric g^. In this case 

m NR (ip) ~ m A(ip) , (2.15) 

leading to M NR ~ ijiq -A(</?oo) and 

Qnr ^ -m ^4'(<Ax>) = -M NR a{ipoo) , (2.16) 



where a((p) := A'{ip)/A{ip) p. Eq. Q2.16 ) agrees with more detailed studies of stellar struc- 



ture in scalar-tensor gravity [25, 26], which also show how Q/M differs from 0(^00) as the 



underlying stellar equation of state becomes more relativistic. 



3 These divergences arise already at the classical level, and can be handled systematically by renormalizing 



the particle action [^7], as is also relevant in other contexts |n| 
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What is the function, m((p), appropriate when the object in question is a nonsingular, 
static black hole? From solutions eqs. ( |2.6D we see that static black hole solutions have a con- 
stant scalar profile, tp = p^, on which the asymptotic metric does not depend. Consequently 
M = m^ifoo) is independent of ifoo — which is consistent with Q = —m'(ipoo) = 0, as found 
above. 

For black holes the conclusion that Tn(<p) must be independent of (p should be very general, 



and should also apply to the coefficients h(ip) and k(ip) of the expansion, eq. (2.11). It is this 
general because it is not necessary to couple to matter at all in order to form a black hole, 
and so the effective action for a black hole should be invariant under the shift symmetry, 



p — > p + constant, of the classical field equations, eqs. fl2.3| ) and (2.4). 
2.3 Jacobson's Miracle Hair-Growth Formula 

For massless scalars the asymptotic value, Pea, of the scalar field can be chosen freely, and is 
not related to the physics of a star or black hole. This is easily understood because far from 
the black hole the geometry becomes flat, and any constant value p = p^ solves the scalar 
field equation, □ ip = 0. However a function linear in time, p = ^t, is another solution to 
□ (p = for flat spacetime, and so one might ask what the black hole solution might be that 
asymptotes to this for large r. 

This question need not be completely academic. For instance, for cosmological applica- 
tions the scalar potential, V(p), need not precisely vanish and might drive a time evolution 
for <p far from the black hole. For small V this evolution could be adiabatic on the time-scales 
relevant to black-hole physics, and so be well-approximated by a slow, approximately linear, 
variation of p at infinity. 

Although the exact solution to this problem is not known, a perturbative solution to 
eqs. ( ^ ) and (|2.4| ) is known with the property that <p — > ip^ + /it as r -> oo |||. This 
solution is given as an expansion about the black-hole solution, 

9^ = 9$ + e 2 9$ + ■■■ , 

p = ^(0) + e(p (l) + e 2^(2) + _ ^ (2 1?) 

where e is a small dimensionless parameter and 

9 ( >^ = ~ (l - ^) dt 2 + TZ ^ M j- r + r^l (2.18) 

is the Schwarzschild metric, p^ = ip^ is a constant (which we can take to vanish with no 
loss of generality) , and 

/ ->r:M 

(2.19) 



t 2GM 



(!) 1 J_l 

^ =2GM ±l0g 

is the leading scalar-field perturbation. Both choices of sign solve eq. (|2.4|), but we will shortly 
see that only p + is physically realistic for our later applications. The back-reaction of the 
metric first arises at C(e 2 ), because equation ( |2.3[ ) is quadratic in dip. 
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Comparing this solution at large r to /it shows that \x = e/2GM and so 

2GM 



ip± = /j,t±2GMfxlog 



1 



+ 0[(GMfj,Y}. (2.20) 



Thus an expansion in e < 1 corresponds to choosing the time scale \x 1 over which the 
asymptotic scalar field ip evolves to be very long relative to the light-crossing time scale, 
2GM/c 3 , of the black hole. 

The necessity for the r-dependence in (p^ 1 ' can be seen by asking how the solution behaves 
near the black-hole event horizon. In ref. (l[ Jacobson shows that the singularity at r = 2GM 
in the spatial part of <pW is required to cancel the singularity due to the breakdown of the 
coordinate t at the horizon. To see this, we re-write it using coordinates that do not break 
down there, such as 

u = t-r*, v = t + r*, (2.21) 



where 



r + 2GM log 



1 



2GM 

is the 'tortoise' coordinate. Trading t for v, the Schwarzschild metric is 



(2.22) 



g$dx»dx u = - ( 1 - ^^-) dv 2 + 2 dr dv + r 2 dQ 2 2 , (2.23) 



and the linearized scalar solutions become 



Now, for an eternal black hole the future event horizon is at t — > oo and r = 2GM, which 
corresponds to u — > oo and finite v. The past horizon is at t — > — oo and r = 2GM, which 
corresponds to v — > — oo and finite u. Therefore tp+ is regular at the future horizon and 
singular at the past horizon, while <p_ is regular at the past horizon and singular at the 
future horizon. Thus, it is the perturbation ip^ that is relevant for a physically realistic 
black hole formed by gravitational collapse with no past horizon. 

Quantifying how much hair 

Because ip+ has a radial nontrivial profile, it also introduces a nonzero scalar charge. Ex- 
panding eq. Q2.19Q in powers of 1/r yields 

and so comparing the 1/r term with —GQ/r implies 

^ = 2e = 4GM/x, (2.26) 
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indicating that asymptotic scalar time-dependence gives black holes nonzero scalar charge. 

Another route to the same conclusion is to recognize that so long as we restrict to spher- 
ically symmetric field configurations with a time- independent metric and dt<p> nonzero but 
t- independent, it follows that U<p = still implies the r-independence of \f—g g r ^d^<p = 
f g<p' s'm9. Using Jacobson's asymptotic expression for large r then shows its value for all 
nonzero r must be 



fir) g(r) <p' = lim f(r) g(r) (p' = lim r V = 2GM e , (2.27) 

r— >oo 



and so GQ = (2GM) 2 fx, as concluded above. What is important for later post-Newtonian 
applications is the black hole's dimensionless charge-to-mass ratio, 

a := j- = 2e = 4GM/J. (2.28) 

We saw above that for ordinary stars with non-relativistic equations of state in quasi Brans- 
Dicke theories this quantity would play the role of the effective coupling, a(</?oo), to the scalar 
field far from the star. This need no longer be true for relativistic systems like neutron stars 
or black holes || 25, ||]. 



The utility of this last derivation lies in its connection to eq. ([2.13 ), which shows that 
the point-particle effective action for such a black hole must satisfy 



-G 

Sip 



f(r)g(r)<p'(r) = (2GM) 2 f,. (2.29) 



Notice that this conclusion is not inconsistent with the black-hole shift symmetry, since this 
symmetry does not preclude a nonzero functional derivative for S p . 

Orbital corrections 

Where the existence of the shift symmetry does play an important role is in how the black 
hole responds to applied scalar and gravitational fields (such as those due to a companion 
within a binary system). For these applications a derivative expansion is appropriate, and 
gives 

S p = J dr \m{ip) + h(<p) z»z v 8 u (p + k((p) dp<p d»(p + • • • ] , (2.30) 

where the functions h(<p) and k(tp) are to be determined by an appropriate matching calcu- 
lation. A possible term linear in zP — i.e. J r dr /(</?) d^tp — is not written in eq. ( p.30| ) 
because it is a total derivative and so does not contribute to the local equations of motion, 
and several other terms involving the background Ricci tensor can be removed by performing 
a judicious field redefinition |2^, [3^]. The ellipses include terms with either more powers of 
d^(f, derivatives of the metric, or higher derivatives of (p. 

As before, for applications to black holes the functions m(ip), h((p) and k(<p) are indepen- 
dent of ip, since black hole solutions do not require couplings between tp and matter and so 
do not break the classical symmetry under constant shifts: <p —> tp + uj. When calculating the 
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response of the black hole to fluctuations, 6<p = p — (foo — [it, about the asymptotic scalar, 
the above action shows the leading terms linear in 5ip are 



8S P ~ 2fi / dr 



~ 2/j J dr 



hi 2 + kg u )5ip tt + hir5p : 



(h-k)5p. t + hf 5p, r + 0(v 2 ) 



(2.31) 



where 5p yT := d r 5p> and 5ip t t '■= fyStp. Notice that the shift symmetry, p — > p + u, guarantees 
that bp always appears differentiated, and so in principle eq. ( 2.3l| ) gives a different response 
than would the expansion 771(92) ~ uiq + m' 5<p + • • • that dominates the response of other 
matter sources, like neutron stars. 

Considerations like these show that computing the changes to black hole orbits due to an 
asymptotic scalar evolution can be complicated, even in the post-Newtonian approximation. 
However the same issues do not similarly complicate the rate that scalars are emitted by 
black holes in binary systems, provided we work to leading nontrivial order in \i. This is 
because the black hole scalar charge is itself of order /j, and so the leading contribution to 
dipole radiation comes from the radiation due to a charge of order /i moving along the zeroth- 
order Schwarzschild orbits. Although the 0(/j 2 ) corrections to these orbits are complicated 
to compute, they are subdominant contributions to the scalar dipole radiation. 

2.4 Dipole radiation from binaries 

Consider first the problem of how a single black hole radiates energy into the scalar field as a 
function of its trajectory, which we assume (temporarily) to be given. Although we assume in 
this section that the scalar is massless, it really suffices that it be light enough to have modes 
into which the system can radiate. Since periodic motion with angular frequency Q typically 



radiates dominantly into modes with angular frequency 0, this requires mc /h <C £1 [33]. 



Scalar radiation from a single particle moving on a given trajectory 

In order to find the power emitted into scalar radiation, the scalar field far from the source 



is written as 



(/j(x,/j) 



G 

lit 

r 



r p y ' t*. 

1* + 7^+2H 



■3 



< + ••• 



! O I -4 



(2.32) 



where V'ii— i; (*) i s the time-dependent STF (symmetric trace- free) multipole moment of order 
I [34|. Primes denote differentiation with respect to time and the right-hand side is to be 
evaluated at the retarded time, t r := t — r, where r := |x| denotes the (large) distance from 
the source to the field point of interest. The terms in the square brackets may be thought of 
as a series in v, the typical velocity of the source. Terms that fall off as 1/r 2 and faster do 
not contribute to the emitted power. The contribution to the emitted power from the dipole 
moment is given by || 

< P 4EWf (2-33) 
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Explicitly, in terms of the non-relativistic source density, g(x,t), defined by the field 
equation 

n f <p = A7rGQ + 0(v 4 ), (2.34) 

where v is the typical velocity of the source, and □/ is the flat-space d'Alembertian operator, 
we have 

ifj(t r ) = J d 3 x g(t r ,x) + 0(v 2 ) 

ipi(t r ) = / d 3 x g(t r ,x) Xi + 0(v 2 ) , (2.35) 



and so on. 

For a point object with scalar charge Q and trajectory z(t) we have g(t, x) = Q <5 3 [x— z(i)], 
and so ip = Q + 0{v 2 ) and ^(t) = Q Zi(t) + 0(v 2 ), etc. Equations Q), and ( p5|) 



show that the Q that appears in standard calculations of the scalar dipole radiation rate is 



the same Q that is identified in section 2.2 using Jacobson's Miracle Hair Growth formula 
for asymptotically varying scalar fields. This allows the use of standard expressions for scalar 
radiation when computing the leading dipole scalar radiation rate for binary systems. 

Scalar radiation from binary systems 

Next consider the radiation rate from a pair of compact objects having masses M A ,M B , and 
scalar charges, Q A , Q B . An important role is played by the ratios a A = Q A /M A , a B = Q B /M B . 
As emphasized earlier, this ratio would be the asymptotic scalar-matter coupling, o((/?oo), 
for ordinary stars in Brans-Dicke-like theories, but this need not be true for objects with 
relativistic structure, like neutron stars or black holes. 

As usual, when all other things are equal, we expect scalar radiation to be well constrained 
by binary orbits because it can arise as dipole radiation. This is by contrast with gravitational 
radiation, which first appears at quadrupole order and so with a rate that is suppressed 
relative to dipole radiation by additional factors of v, with v being a representative orbital 
speed (which for numerical purposes we take to be the average orbital speed). 

To quantify this, let F mon , F dip and F quad respectively denote the power in monopole, 
dipole and quadrupole radiation, integrated over an orbital period. For objects whose orbits 
are close to Newtonian the total energy loss is the sum of tensor and scalar contributions, 
F = F g + F^, with the leading gravitational contribution being quadrupole radiation, F g ~ 
^quad _|_..._ By contrast, for scalars monopole and dipole radiation is possible, F v ~ F™ OTi + 

P^V£> _|_ ^ uad _|_ . . . 

In order of magnitude F™ n = 0[a\ B v% Ff v = 0[a 2 A B v% i^ quad = 0[v 5 } and F£ uad = 
0[a 2 A B v 5 ]. (Notice that F™ on is order v 5 rather than order v because the scalar charges Q AiB 
are time- independent.) Constraints on scalar-tensor theories arise as upper bounds on the 
scalar radiation rate for systems where the gravitational rate is observed. To this end define 
the ratio 



7? 7? di P 



S := ye. ~ ±* = o 

F g F^ nad 



v 2 



(2.36) 
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Explicit formulae for the orbitally averaged power emitted into scalar and gravity waves 
are calculated in chapter 6 of ||, which finds (re- introducing the factors of c) 

= y,,fiM ^^- )! (i + «]) , (,ar, 

*r - 32 ^ (CM " ) ^;^ +37e4/a6) ( 1 + «i) . (-38) 

where e is the orbit eccentricity, Q = 2tt/P where P is the orbital period, M = M A + M B is 
the total mass, and v = M A M B /M 2 . Using these, the ratio £ becomes 

_ 5c 2 (a A -a a ) 2 (l-e 2 )(l + e 2 /2) 
96(GMft) 2 / 3 (l + 73e 2 /24 + 37e 4 /96) 
5 (a A - a fl ) 2 



1 + 0[(v/c) 

-96 (,/c) 2 /W^ + ^/cnj. (2-39) 

Notice that (GMS7/c 3 ) 2 / 3 = GM/ac 2 = (v/c) 2 in the Keplerian limit, where o is the semi- 
major axis of the orbit and v is the average orbital speed. 

For two black holes with an asymptotically static scalar field we have a A = a B = because 
of the no-hair condition Q A = Q B = 0, as is well-known to hold for Kerr and Schwarzschild 
black holes in general relativity [p5|] . This implies in particular the usual absence of dipole 
radiation, £ ~ 0, as expected. 

When the asymptotic scalar varies slowly in time — i.e. fi ^ — we've seen that the 
above expressions can also be used to compute the leading contribution to the scalar dipole 



radiation, with the scalar charge given as in section |2.3| : a A = 2e A and a B = 2e B . Since 
both black holes see the same asymptotic scalar field far from the binary system, consistency 
requires 

wk=^k = "' <2 - 40) 

and so because the two black holes share the same asymptotic time-dependent scalar field, 
their scalar charges are related to one another. 

The quantity appearing in £ is a A — a B = 2(e A — e B ) = 4Gfi(M A — M B ). We see that 
inspiralling black holes with differing masses will radiate dipole radiation when situated within 
a slowly varying asymptotic scalar field, provided only that the scalar is light enough to 
receive the radiation. The amount of radiation is guaranteed to be small by virtue of the 
approximation e A)B -C 1 that underlies our linearized analysis. 

This calculation assumes Q A and M A are fixed quantities over any one orbit, and this is 
true in the present case if e AjB <C 1 and [i doesn't change in time. Notice in particular that 
(at leading order) we need not assume [i <C £1 for the radiation-loss calculation to be valid. 

3. The Quasar OJ287 



We next apply the formalism of section 2.4 to OJ287. This is a quasar situated at redshift 



z = 0.306 that undergoes periodic bursts. This system has been argued to consist of two 
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Parameter 


Meaning 


Value 




Precession rate 


39.1(1)° per orbit 


M A 


Mass of black hole A 


1.84(1) • 1O 1O M 


M B 


Mass of black hole B 


1.40(3) • 1O 8 M 


Xa 


Dimensionless Kerr parameter of black hole A 


0.28(3) 


e 


Orbital eccentricity 


0.658(1) 


P 


Orbital period 


~ 9 years ~ 2.8 • 10 8 s 


a 


Semi- major axis 


~ 1.7- 10 15 m 


^orb 


Typical orbital velocity 


~ 0.13c 



Table 1: Orbital parameters of OJ287. These are the 'intrinsic' values, as measured in the centre-of- 
mass of the binary. The values measured on Earth differ by a redshift factor. Data taken from p0[ . 
Uncertainties are given at the 3 sigma level. 



gravitationally-bound black holes [36, 37, 38], with the bursts occurring when the 'small' 
black hole (black hole B, with mass M B ~ 10 8 M©) passes through the accretion disc of the 
larger one (black hole A, which is a monster with M A ~ 10 10 Mq). 

The orbital parameters of the system have been reconstructed from the observed bursting 
pattern, leading to the orbital parameters summarized in table In particular, the success of 
the orbital description requires the orbital decay due to the emission of gravitational radiation, 
with the observed radiated power being within 6% of the prediction of General Relativity [37]. 

Taking this binary black hole interpretation as correct, we may use its success to infer 
an upper limit for the time-variation of any ambient light scalar fields that were present in 
the black hole neighborhood at an epoch z = 0.306. Notice that this does not require making 
any assumptions about the nature of the scalar-matter couplings. This bound applies to any 
scalar fields at this epoch, provided that they are light enough that it is possible to radiate 
into their scalar modes. This requires the mass of the scalar to be much lighter than the 
angular frequency, f2. Given that the orbital period is about 9 years, the corresponding scalar 
mass limit becomes m >C 10~ 23 eV/c 2 [33]. 

Using the black hole orbital parameters in equation ( 12.39 ) , we find £ ~ 0.94(a A — a B ) 2 , 
and so imposing £ < 0.06 gives the constraint \a A — a B \ < 0.25 and so \e A — e B \ < 0.13. We 
have e AjB = 2GM A>B fi/c s , where 2GM A /c 3 ~ 1.8 • 10 5 s ~ 2 day is the light -crossing time of 
black hole A, and 2GM B /c 3 ~ 1.4 • 10 3 s ~ 23 min is that of black hole B. 

This finally implies > 1.4 • 10 6 s ~ 16 days. This limit holds even though it is 

much smaller than the orbital period, because the instantaneous radiation rate relies only on 
the constancy of [i and the black hole masses and scalar charges. Our upper bound on 
implies |e A | < 0.13 and |e s | < 1.0 • 10~ 3 . This is consistent with our linearization assumption 
e A , B < 1. 



Notice also that the spin of black hole A is claimed to have been measured [39]. The 
solution described in section 2.3 can be generalized to rotating black holes M, and the analogue 
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of equation fl2.28|) for Kerr black holes is 

a A = 2GM A fi(l + y/l-xl) , (3-1) 

where \a = cJ A jGM\ is the dimensionless Kerr parameter, and J A is the angular momentum 
of the black hole. It follows from equation (|3.l| ) that taking spin into account does not 
significantly change our bounds. 

4. Summary and Outlook 

In this paper we use Jacobson's Miracle Hair-Growth Formula to show that inspiralling black 
holes can radiate dipole scalar radiation whenever there exists a scalar whose mass is small 
enough to be radiated (typically this requires the mass to be much smaller than the orbital 
period, m <C £2). We then apply this radiation calculation to the quasar OJ287, which has 
been interpreted as a very massive double black hole binary. Using the inferred black hole 
masses and orbital properties, and the success of the GR prediction for the orbital decay rate, 
we conclude that any ambient scalar fields lighter than m < 10 -23 eV at this epoch (redshift 
z = 0.306) must have instantaneous time evolution over scales larger than •> 16 days. 

Although this is not a spectacular bound compared with cosmological time evolution, 
there are two things about it that are quite remarkable. First, it is completely independent of 
how the light scalar couples to matter since this does not enter into the amount of acquired 
black hole hair. Second, it directly constrains the instantaneous scalar evolution rate, rather 
than an average change in scalar field values over cosmological time intervals. 

This system is representative of the rich kinds of physics that can arise for gravity/scalar 
systems with time-dependent environments, and indicates the wealth of information that is 
likely to become available once gravitational radiation is detected from astrophysical sources. 

Direct detection of gravitational waves in the future will provide a systematic means of 
measuring the properties of double black hole binaries, making it possible to map out the 
bounds on at different redshifts. In preparation for these happy times it would be 

useful to predict in more detail the response of these systems to scalar-tensor theories, and 
in particular how their orbits respond to their scalar 'environment'. 

Alternatively, if the scalars are not precisely massless they will propagate with different 
speeds than those of the primary gravitational signal, allowing potential time delays in grav- 
itational wave signals whose presence could signal the existence of a new scalar component 
to gravity. 
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